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1. 


INTRODUCTION 


The  Summer  Research  Program  (SRP),  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
(AFOSR),  offers  paid  opportunities  for  university  faculty,  graduate  students,  and  high  school 
students  to  conduct  research  in  U.S.  Air  Force  research  laboratories  nationwide  during  the  summer. 

Introduced  by  AFOSR  in  1978,  this  innovative  program  is  based  on  the  concept  of  teaming 
academic  researchers  with  Air  Force  scientists  in  the  same  disciplines  using  laboratory  facilities  and 
equipment  not  often  available  at  associates'  institutions. 

The  Summer  Faculty  Research  Program  (SFRP)  is  open  annually  to  approximately  150  faculty 
members  with  at  least  two  years  of  teaching  and/or  research  experience  in  accredited  U.S.  colleges, 
universities,  or  technical  institutions.  SFRP  associates  must  be  either  U.S.  citizens  or  permanent 
residents. 

The  Graduate  Student  Research  Program  (GSRP)  is  open  annually  to  approximately  100  graduate 
students  holding  a  bachelor's  or  a  master's  degree;  GSRP  associates  must  be  U.S.  citizens  enrolled 
full  time  at  an  accredited  institution. 

The  High  School  Apprentice  Program  (HSAP)  annually  selects  about  125  high  school  students 
located  within  a  twenty  mile  commuting  distance  of  participating  Air  Force  laboratories. 

AFOSR  also  offers  its  research  associates  an  opportunity,  under  the  Summer  Research  Extension 
Program  (SREP),  to  continue  their  AFOSR-sponsored  research  at  their  home  institutions  through  the 
award  of  research  grants.  In  1994  the  maximum  amount  of  each  grant  was  increased  from  $20,000 
to  $25,000,  and  the  number  of  AFOSR-sponsored  grants  decreased  from  75  to  60.  A  separate 
annual  report  is  compiled  on  the  SREP. 

The  numbers  of  projected  summer  research  participants  in  each  of  the  three  categories  and  SREP 
“grants”  are  usually  increased  through  direct  sponsorship  by  participating  laboratories. 

AFOSR's  SRP  has  well  served  its  objectives  of  building  critical  links  between  Air  Force  research 
laboratories  and  the  academic  community,  opening  avenues  of  communications  and  forging  new 
research  relationships  between  Air  Force  and  academic  technical  experts  in  areas  of  national 
interest,  and  strengthening  the  nation's  efforts  to  sustain  careers  in  science  and  engineering.  The 
success  of  the  SRP  can  be  gauged  from  its  growth  from  inception  (see  Table  1)  and  from  the 
favorable  responses  the  1997  participants  expressed  in  end-of-tour  SRP  evaluations  (Appendix  B). 

AFOSR  contracts  for  administration  of  the  SRP  by  civilian  contractors.  The  contract  was  first 
awarded  to  Research  &  Development  Laboratories  (RDL)  in  September  1990.  After  completion  of 
the  1990  contract,  RDL  (in  1993)  won  the  recompetition  for  the  basic  year  and  four  1-year  options. 


1 


2. 


PARTICIPATION  IN  THE  SUMMER  RESEARCH  PROGRAM 


The  SRP  began  with  faculty  associates  in  1979;  graduate  students  were  added  in  1982  and  high 
school  students  in  1986.  The  following  table  shows  the  number  of  associates  in  the  program  each 
year. 


Beginning  in  1993,  due  to  budget  cuts,  some  of  the  laboratories  weren’t  able  to  afford  to  fund  as 
many  associates  as  in  previous  years.  Since  then,  the  number  of  funded  positions  has  remained 
fairly  constant  at  a  slightly  lower  level. 


3.  RECRUITING  AND  SELECTION 

The  SRP  is  conducted  on  a  nationally  advertised  and  competitive-selection  basis.  The  advertising 
for  faculty  and  graduate  students  consisted  primarily  of  the  mailing  of  8,000  52-page  SRP  brochures 
to  chairpersons  of  departments  relevant  to  AFOSR  research  and  to  administrators  of  grants  in 
accredited  universities,  colleges,  and  technical  institutions.  Historically  Black  Colleges  and 
Universities  (HBCUs)  and  Minority  Institutions  (Mis)  were  included.  Brochures  also  went  to  all 
participating  USAF  laboratories,  the  previous  year’s  participants,  and  numerous  individual 
requesters  (over  1000  annually). 

RDL  placed  advertisements  in  the  following  publications:  Black  Issues  in  Higher  Education,  Winds 
of  Change,  and  IEEE  Spectrum.  Because  no  participants  list  either  Physics  Today  or  Chemical  & 
Engineering  News  as  being  their  source  of  learning  about  the  program  for  the  past  several  years, 
advertisements  in  these  magazines  were  dropped,  and  the  funds  were  used  to  cover  increases  in 
brochure  printing  costs. 

High  school  applicants  can  participate  only  in  laboratories  located  no  more  than  20  miles  from  their 
residence.  Tailored  brochures  on  the  HSAP  were  sent  to  the  head  counselors  of  180  high  schools  in 
the  vicinity  of  participating  laboratories,  with  instructions  for  publicizing  the  program  in  their 
schools.  High  school  students  selected  to  serve  at  Wright  Laboratory's  Armament  Directorate 
(Eglin  Air  Force  Base,  Florida)  serve  eleven  weeks  as  opposed  to  the  eight  weeks  normally  worked 
by  high  school  students  at  all  other  participating  laboratories. 

Each  SFRP  or  GSRP  applicant  is  given  a  first,  second,  and  third  choice  of  laboratory.  High  school 
students  who  have  more  than  one  laboratory  or  directorate  near  their  homes  are  also  given  first, 
second,  and  third  choices. 

Laboratories  make  their  selections  and  prioritize  their  nominees.  AFOSR  then  determines  the 
number  to  be  funded  at  each  laboratory  and  approves  laboratories'  selections. 

Subsequently,  laboratories  use  their  own  funds  to  sponsor  additional  candidates.  Some  selectees  do 
not  accept  the  appointment,  so  alternate  candidates  are  chosen.  This  multi-step  selection  procedure 
results  in  some  candidates  being  notified  of  their  acceptance  after  scheduled  deadlines.  The  total 
applicants  and  participants  for  1998  are  shown  in  this  table. 
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1998  Applicants  and  Participants 

PARTICIPANT 

TOTAL 

SELECTEES 

DECLINING 

CATEGORY 

APPLICANTS 

SELECTEES 

SFRP 

382 

85 

13 

(HBCU/MI) 

(0) 

(0) 

(0) 

GSRP 

130 

40 

7 

(HBCU/MI) 

(0) 

(0) 

(0) 

HSAP 

328 

88 

22 

TOTAL 

840 

213 

42 

4.  SITE  VISITS 

During  June  and  July  of  1998,  representatives  of  both  AFOSR/NI  and  RDL  visited  each 
participating  laboratory  to  provide  briefings,  answer  questions,  and  resolve  problems  for  both 
laboratory  personnel  and  participants.  The  objective  was  to  ensure  that  the  SRP  would  be  as 
constructive  as  possible  for  all  participants.  Both  SRP  participants  and  RDL  representatives  found 
these  visits  beneficial.  At  many  of  the  laboratories,  this  was  the  only  opportunity  for  all  participants 
to  meet  at  one  time  to  share  their  experiences  and  exchange  ideas. 


5.  HISTORICALLY  BLACK  COLLEGES  AND  UNIVERSITIES  AND  MINORITY 
INSTITUTIONS  (HBCU/MIs) 

Before  1993,  an  RDL  program  representative  visited  from  seven  to  ten  different  HBCU/MIs 
annually  to  promote  interest  in  the  SRP  among  the  faculty  and  graduate  students.  These  efforts  were 
marginally  effective,  yielding  a  doubling  of  HBCI/MI  applicants.  In  an  effort  to  achieve  AFOSR’s 
goal  of  10%  of  all  applicants  and  selectees  being  HBCU/MI  qualified,  the  RDL  team  decided  to  try 
other  avenues  of  approach  to  increase  the  number  of  qualified  applicants.  Through  the  combined 
efforts  of  the  AFOSR  Program  Office  at  Bolling  AFB  and  RDL,  two  very  active  minority  groups 
were  found,  HACU  (Hispanic  American  Colleges  and  Universities)  and  AISES  (American  Indian 
Science  and  Engineering  Society).  RDL  is  in  communication  with  representatives  of  each  of  these 
organizations  on  a  monthly  basis  to  keep  up  with  the  their  activities  and  special  events.  Both 
organizations  have  widely-distributed  magazines/quarterlies  in  which  RDL  placed  ads. 

Since  1994  the  number  of  both  SFRP  and  GSRP  HBCU/MI  applicants  and  participants  has  increased 
ten-fold,  from  about  two  dozen  SFRP  applicants  and  a  half  dozen  selectees  to  over  100  applicants 
and  two  dozen  selectees,  and  a  half-dozen  GSRP  applicants  and  two  or  three  selectees  to  18 
applicants  and  7  or  8  selectees.  Since  1993.  the  SFRP  had  a  two-fold  applicant  increase  and  a  two¬ 
fold  selectee  increase.  Since  1993.  the  GSRP  had  a  three-fold  applicant  increase  and  a  three  to  four¬ 
fold  increase  in  selectees. 
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In  addition  to  RDL's  special  recruiting  efforts,  APOSR  attempts  each  year  to  obtain  additional 
funding  or  use  leftover  funding  from  cancellations  the  past  year  to  fund  HBCU/MI  associates. 


SRP  HBCU/MI  Participation,  By  Year 

YEAR 

SFRP 

GSRP 

Applicants 

Participants 

Applicants 

Participants 

1985 

76 

23 

15 

11 

1986 

70 

18 

20 

10 

1987 

82 

32 

32 

10 

1988 

53 

17 

23 

14 

1989 

39 

15 

13 

4 

1990 

43 

14 

17 

3 

1991 

42 

13 

8 

5 

1992 

70 

13 

9 

5 

1993 

60 

13 

6 

2 

1994 

90 

16 

11 

6 

1995 

90 

21 

20 

8 

1996 

119 

27 

18 

7 

6.  SRP  FUNDING  SOURCES 

Funding  sources  for  the  1998  SRP  were  the  AFOSR-provided  slots  for  the  basic  contract  and 
laboratory  funds.  Funding  sources  by  category  for  the  1998  SRP  selected  participants  are  shown 
here. 
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1998  SRP  FUNDING  CATEGORY 

SFRP 

GSRP 

HSAP 

AFOSR  Basic  Allocation  Funds 

67 

38 

75 

USAF  Laboratory  Funds 

17 

2 

13 

Slots  Added  by  AFOSR 

(Leftover  Funds) 

0 

0 

0 

HBCU/MI  By  AFOSR 
(Using  Procured  Addn’l  Funds) 

0 

0 

N/A 

TOTAL 

84 

40 

88 

7.  COMPENSATION  FOR  PARTICIPANTS 

Compensation  for  SRP  participants,  per  five-day  work  week,  is  shown  in  this  table. 


1998  SRP  Associate  Compensation 


PARTICIPANT  CATEGORY 

1991 

1992 

1993 

1994 

1995 

1996 

1997 

1998 

Faculty  Members 

$690 

$718 

$740 

$740 

$740 

$770 

$770 

$793 

Graduate  Student 
(Master's  Degree) 

$425 

$442 

$455 

$455 

$455 

$470 

$470 

$484 

Graduate  Student 

(Bachelor's  Degree) 

$365 

$380 

$391 

$391 

$391 

$400 

$400 

$412 

High  School  Student 
(First  Year) 

$200 

$200 

$200 

$200 

$200 

$200 

$200 

$200 

High  School  Student 
(Subsequent  Years) 

$240 

$240 

$240 

$240 

$240 

$240 

$240 

$240 

The  program  also  offered  associates  whose  homes  were  more  than  50  miles  from  the  laboratory  an 
expense  allowance  (seven  days  per  week)  of  $52/day  for  faculty  and  $41 /day  for  graduate  students. 
Transportation  to  the  laboratory  at  the  beginning  of  their  tour  and  back  to  their  home  destinations  at 
the  end  was  also  reimbursed  for  these  participants.  Of  the  combined  SFRP  and  GSRP  associates, 

65  °Ic  claimed  travel  reimbursements  at  an  average  round-trip  cost  of  $730. 

Faculty  members  were  encouraged  to  visit  their  laboratories  before  their  summer  tour  began.  All 
costs  of  these  orientation  visits  were  reimbursed.  Forty-three  percent  (85  out  of  188)  of  faculty 
associates  took  orientation  trips  at  an  average  cost  of  $449.  By  contrast,  in  1993,  58  %  of  SFRP 
associates  elected  to  take  an  orientation  visits  at  an  average  cost  of  $685;  that  was  the  highest 
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percentage  of  associates  opting  to  take  an  orientation  trip  since  RDL  has  administered  the  SRP,  and 
the  highest  average  cost  of  an  orientation  trip. 

Program  participants  submitted  biweekly  vouchers  countersigned  by  their  laboratory  research  focal 
point,  and  RDL  issued  paychecks  so  as  to  arrive  in  associates'  hands  two  weeks  later. 

This  is  the  third  year  of  using  direct  deposit  for  the  SFRP  and  GSRP  associates.  The  process  went 
much  more  smoothly  with  respect  to  obtaining  required  information  from  the  associates,  about  15% 
of  the  associates’  information  needed  clarification  in  order  for  direct  deposit  to  properly  function  as 
opposed  to  7%  from  last  year.  The  remaining  associates  received  their  stipend  and  expense 
payments  via  checks  sent  in  the  US  mail. 

HSAP  program  participants  were  considered  actual  RDL  employees,  and  their  respective  state  and 
federal  income  tax  and  Social  Security  were  withheld  from  their  paychecks.  By  the  nature  of  their 
independent  research,  SFRP  and  GSRP  program  participants  were  considered  to  be  consultants  or 
independent  contractors.  As  such,  SFRP  and  GSRP  associates  were  responsible  for  their  own 
income  taxes,  Social  Security,  and  insurance. 

8.  CONTENTS  OF  THE  1998  REPORT 

The  complete  set  of  reports  for  the  1998  SRP  includes  this  program  management  report  (Volume  1) 
augmented  by  fifteen  volumes  of  final  research  reports  by  the  1998  associates,  as  indicated  below: 


1998  SRP  Final  Report  Volume  Assignments 


LABORATORY 

SFRP 

GSRP 

HSAP 

Armstrong 

2 

7 

12 

Phillips 

3 

8 

13 

Rome 

4 

9 

14 

Wright 

5A,  5B 

10 

15 

AEDC,  ALCs,  USAFA,  WHMC 

6 

11 

7 


APPENDIX  A  -  PROGRAM  STATISTICAL  SUMMARY 


A.  Colleges/Universities  Represented 

Selected  SFRP  associates  represented  169  different  colleges,  universities,  and  institutions, 
GSRP  associates  represented  95  different  colleges,  universities,  and  institutions. 

B.  States  Represented 

SFRP  -Applicants  came  from  47  states  plus  Washington  D.C.  Selectees  represent  44  states. 
GSRP  -  Applicants  came  from  44  states.  Selectees  represent  32  states. 

HSAP  -  Applicants  came  from  thirteen  states.  Selectees  represent  nine  states. 


Total  Number  of  Participants 

SFRP 

85 

GSRP 

40 

HSAP 

88 

TOTAL 

213 

Degrees  Represented 

SFRP 

GSRP 

TOTAL 

Doctoral 

83 

0 

83 

Master's 

1 

3 

4 

Bachelor's 

0 

22 

22 

TOTAL 

186 

25 

109 

A-l 


SFRP  Academic  Titles 

Assistant  Professor 

36 

Associate  Professor 

34 

Professor 

15 

Instructor 

0 

Chairman 

0 

Visiting  Professor 

0 

Visiting  Assoc.  Prof. 

0 

Research  Associate 

0 

TOTAL 

85 

Source  of  Learning  About  the  SRP 


Category 

Applicants 

Selectees 

Applied/participated  in  prior  years 

177 

47 

Colleague  familiar  with  SRP 

104 

24 

Brochure  mailed  to  institution 

101 

21 

Contact  with  Air  Force  laboratory 

101 

39 

IEEE  Spectrum 

12 

1 

BIIHE 

4 

0 

Other  source 

117 

30 

TOTAL 


616 


162 


APPENDIX  B  -  SRP  EVALUATION  RESPONSES 


1.  OVERVIEW 

Evaluations  were  completed  and  returned  to  RDL  by  four  groups  at  the  completion  of  the  SRP.  The 
number  of  respondents  in  each  group  is  shown  below. 


Table  B-l .  Total  SRP  Evaluations  Received 


Evaluation  Group 

Responses 

SFRP  &  GSRPs 

100 

HSAPs 

75 

USAF  Laboratory  Focal  Points 

84 

USAF  Laboratory  HSAP  Mentors 

6 

All  groups  indicate  unanimous  enthusiasm  for  the  SRP  experience. 


The  summarized  recommendations  for  program  improvement  from  both  associates  and  laboratory 
personnel  are  listed  below: 

A.  Better  preparation  on  the  labs’  pan  prior  to  associates'  arrival  (i.e.,  office  space, 
computer  assets,  clearly  defined  scope  of  work). 

B.  Faculty  Associates  suggest  higher  stipends  for  SFRP  associates. 

C.  Both  HSAP  Air  Force  laboratory  mentors  and  associates  would  like  the  summer  tour 

extended  from  the  current  8  weeks  to  either  10  or  11  weeks;  the  groups  state  it  takes 
4-6  weeks  just  to  get  high  school  students  up-to-speed  on  what’s  going  on  at 

laboratory.  (Note:  this  same  argument  was  used  to  raise  the  faculty  and  graduate 

student  participation  time  a  few  years  ago.) 
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2.  1998  USAF  LABORATORY  FOCAL  POINT  (LFP)  EVALUATION  RESPONSES 


The  summarized  results  listed  below  are  from  the  84  LFP  evaluations  received. 
1.  LFP  evaluations  received  and  associate  preferences: 

Table  B-2.  Air  Force  LFP  Evaluation  Responses  (By  Type) 


How  Many  Associates  Would  You  Prefer  To  Get  ? _ (%  Response) 


Lab 

EvaLs 

Reev’d 

0 

SFRP 

1  2 

3+ 

GSRP  (w/Univ  Professor) 

0  12  3+ 

GSRP  (w/o  Univ  Professor) 

0  12  3+ 

AEDC 

0 

- 

- 

- 

- 

- 

- 

- 

- 

- 

_ 

_ 

.. 

WHMC 

0 

- 

- 

- 

- 

- 

- 

- 

- 

_ 

- 

- 

- 

AL 

7 

28 

28 

28 

14 

54 

14 

28 

0 

86 

0 

14 

0 

USAFA 

1 

0 

100 

0 

0 

100 

0 

0 

0 

0 

100 

0 

0 

PL 

25 

40 

40 

16 

4 

88 

12 

0 

0 

84 

12 

4 

0 

RL 

5 

60 

40 

0 

0 

80 

10 

0 

0 

100 

0 

0 

0 

WL 

46 

30 

43 

20 

6 

78 

17 

4 

0 

93 

4 

2 

0 

Total 

84 

32% 

50% 

13% 

5% 

80% 

11% 

6% 

0% 

73% 

23% 

4% 

0% 

LFP  Evaluation  Summary.  The  summarized  responses,  by  laboratory,  are  listed  on  the  following 
page.  LFPs  were  asked  to  rate  the  following  questions  on  a  scale  from  1  (below  average)  to  5 
(above  average). 

2.  LFPs  involved  in  SRP  associate  application  evaluation  process: 

a.  Time  available  for  evaluation  of  applications: 

b.  Adequacy  of  applications  for  selection  process: 

3.  Value  of  orientation  trips: 

4.  Length  of  research  tour: 

5  a.  Benefits  of  associate’s  work  to  laboratory: 
b.  Benefits  of  associate's  work  to  Air  Force: 

6.  a.  Enhancement  of  research  qualifications  for  LFP  and  staff: 

b.  Enhancement  of  research  qualifications  for  SFRP  associate: 

c.  Enhancement  of  research  qualifications  for  GSRP  associate: 

7.  a.  Enhancement  of  knowledge  for  LFP  and  staff: 

b.  Enhancement  of  knowledge  for  SFRP  associate: 

c.  Enhancement  of  knowledge  for  GSRP  associate: 

8.  Value  of  Air  Force  and  university  links: 

9.  Potential  for  future  collaboration: 

10.  a.  Your  working  relationship  with  SFRP: 
b.  Your  working  relationship  w  ith  GSRP: 

11.  Expenditure  of  your  time  worthwhile: 

(Continued  on  next  page) 


B-2 


12.  Quality  of  program  literature  for  associate: 

13.  a.  Quality  of  RDL’s  communications  with  you: 

b.  Quality  of  RDL's  communications  with  associates: 

14.  Overall  assessment  of  SRP: 


Table  B-3.  Laboratory  Focal  Point  Reponses  to  above  questions 


AEDC 

AL 

USAF 

A 

PL 

RL 

WHMC 

WL 

ft  Evals  Reev’d 

7 

1 

14 

5 

0 

46 

Question  # 

2 

- 

86  % 

80  % 

- 

85  % 

2a 

- 

4.3 

4.0 

- 

3.6 

2b 

- 

4.0 

4.5 

- 

4.1 

3 

- 

4.5 

4.3 

- 

3.7 

4 

- 

4.1 

4.2 

- 

3.9 

5a 

- 

4.3 

4.6 

- 

4.4 

5b 

- 

4.5 

n/a 

4.2 

4.6 

- 

4.3 

6a 

- 

4.5 

5.0 

4.0 

4.4 

- 

4.3 

6b 

- 

4.3 

n/a 

4.1 

5.0 

- 

4.4 

6c 

- 

3.7 

5.0 

3.5 

5.0 

- 

4.3 

7a 

- 

4.7 

5.0 

4.0 

4.4 

- 

4.3 

7b 

- 

4.3 

n/a 

4.2 

5.0 

- 

4.4 

7c 

- 

4.0 

5.0 

3.9 

5.0 

- 

4.3 

8 

- 

4.6 

4.0 

4.5 

4.6 

- 

4.3 

9 

- 

4.9 

5.0 

4.4 

4.8 

- 

4.2 

10a 

- 

5.0 

n/a 

4.6 

4.6 

- 

4.6 

10b 

- 

4.7 

5.0 

3.9 

5.0 

- 

4.4 

11 

- 

4.6 

5.0 

4.4 

4.8 

- 

4.4 

12 

- 

4.0 

4.0 

4.0 

4.2 

- 

3.8 

13a 

- 

3.2 

4.0 

3.5 

3.8 

- 

3.4 

13b 

- 

3.4 

4.0 

3.6 

4.5 

- 

3.6 

14 

- 

4.4 

5.0 

4.4 

4.8 

- 

4.4 
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3.  1998  SFRP  &  GSRP  EVALUATION  RESPONSES 


The  summarized  results  listed  below  are  from  the  120  SFRP/GSRP  evaluations  received. 

Associates  were  asked  to  rate  the  following  questions  on  a  scale  from  1  (below  average)  to  5  (above 
average)  -  by  Air  Force  base  results  and  over-all  results  of  the  1998  evaluations  are  listed  after  the 
questions. 

1.  The  match  between  the  laboratories  research  and  your  field: 

2.  Your  working  relationship  with  your  LFP: 

3.  Enhancement  of  your  academic  qualifications: 

4.  Enhancement  of  your  research  qualifications: 

5.  Lab  readiness  for  you:  LFP,  task,  plan: 

6.  Lab  readiness  for  you:  equipment,  supplies,  facilities: 

7.  Lab  resources: 

8.  Lab  research  and  administrative  support: 

9.  Adequacy  of  brochure  and  associate  handbook: 

10.  RDL  communications  with  you: 

11.  Overall  payment  procedures: 

12.  Overall  assessment  of  the  SRP: 

13.  a.  Would  you  apply  again? 

b.  Will  you  continue  this  or  related  research? 

14.  Was  length  of  your  tour  satisfactory? 

15.  Percentage  of  associates  who  experienced  difficulties  in  finding  housing: 

16.  Where  did  you  stay  during  your  SRP  tour? 


a. 

At  Home: 

b. 

With  Friend: 

c. 

On  Local  Economy: 

d. 

Base  Quarters: 

Value 

of  orientation  visit: 

a. 

Essential: 

b. 

Convenient: 

c. 

Not  Worth  Cost: 

d. 

Not  Used: 

SFRP  and  GSRP  associate’s  responses  are  listed  in  tabular  format  on  the  following  page. 
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Table  B-4.  1997  SFRP  &  GSRP  Associate  Responses  to  SRP  Evaluation 


Arnold 

Brooks 

Eglin 

Griffis 

Hanscom 

Kelly 

Rutland 

Lackland 

Robins 

Tyndall 

WPAFB 

averag 

e 

# 

res 

■ 

48 

6 

14 

31 

19 

3 

32 

i 

2 

10 

85 

257 

i 

4.8 

eh 

■Hi 

ESI 

eh 

4.9 

HI 

mm 

5.0 

4.0 

EH 

2 

5.0 

EH 

4.1 

EH 

EH 

4.7 

5.0 

mm 

5.0 

Efikl 

4.8  ! 

ESI 

3 

EH 

4.0 

gin 

4.2 

HI 

mm 

5.0 

5.0 

itfn 

■ 

EH 

4 

3.8 

gfi 

4.4 

5.0 

4.0 

EH 

■Hi 

KEK 

5 

I 

3.3 

4.8 

eh 

4.5 

EH 

mm 

5.0 

3.9 

WEm. 

EH 

6 

Kill 

■SE 

ESI 

EH 

4.5 

3.8 

5.0 
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4.  1998  USAF  LABORATORY  HSAP  MENTOR  EVALUATION  RESPONSES 

Not  enough  evaluations  received  (5  total)  from  Mentors  to  do  useful  summary. 
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5.  1998  HSAP  EVALUATION  RESPONSES 


The  summarized  results  listed  below  are  from  the  23  HSAP  evaluations  received. 

HSAP  apprentices  were  asked  to  rate  the  following  questions  on  a  scale  from 
1  (below  average)  to  5  (above  average) 

1.  Your  influence  on  selection  of  topic/type  of  work. 

2.  Working  relationship  with  mentor,  other  lab  scientists. 

3.  Enhancement  of  your  academic  qualifications. 

4.  Technically  challenging  work. 

5.  Lab  readiness  for  you:  mentor,  task,  work  plan,  equipment. 

6.  Influence  on  your  career. 

7.  Increased  interest  in  math/science. 

8.  Lab  research  &  administrative  support. 

9.  Adequacy  of  RDL’s  Apprentice  Handbook  and  administrative  materials. 

10.  Responsiveness  of  RDL  communications. 

11.  Overall  payment  procedures. 

12.  Overall  assessment  of  SRP  value  to  you. 

13.  Would  you  apply  again  next  year?  Yes  (92  %) 

14.  Will  you  pursue  future  studies  related  to  this  research?  Yes  (68  %) 

15.  Was  Tour  length  satisfactory?  Yes  (82  %) 
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WIND  VALIDATION: 

INCOMPRESSIBLE  TURBULENT  FLOW  PAST  A  FLAT  PLATE 


Gregory  Michael  Laskowski 
Graduate  Student 

Department  of  Aeronautics  and  Astronautics 
Stanford  University 


Abstract 


The  NPARC  alliance  has  invested  considerable  resources  in  developing  the  WIND  code:  a  robust, 
user  friendly,  finite-volume,  structured,  multi-zone,  compressible  flow  solver  with  flexible  chemistry  and 
turbulence  models  which  can  be  run  on  a  variety  of  platforms.  A  validation  effort  was  undertaken  to 
examine  the  ability  of  the  WIND  code’s  accuracy  and  efficiency  for  incompressible  turbulent  flat  plate 
flow.  The  validation  work  presented  herein  presents  the  results  obtained  upon  conducting  a  grid 
independence  study  with  and  without  wall  functions  and  validating  the  different  turbulence  models, 
explicit  Euler  flux  schemes,  order  of  the  Euler  flux  schemes,  implicit  methods,  and  time  integration 
schemes.  Also,  results  obtained  in  conducting  2D  simulations  were  than  compared  with  those  obtained 
for  3D  simulations.  Overall,  the  different  methods  proved  quite  successful,  with  only  slight  problems 
observed  with  some  of  the  test  configurations. 
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WIND  VALIDATION: 

INCOMPRESSIBLE  TURBULENT  FLOW  PAST  A  FLAT  PLATE 

Gregory  Michael  Laskowski 

Introduction 

It  is  generally  believed  in  the  CFD  community  that  the  first  major  breakthrough  in  CFD  can  be 
attributed  to  Murmon  and  Cole’s  solution  of  the  Transonic  Small  Disturbance  equations’,  which  govern 
steady,  irrotational,  isentropic  flow.  Twenty-five  years  later,  through  the  collaborative  and  individual 
efforts  of  academia,  industry  and  the  military,  it  is  now  possible  to  solve  the  full  Navier-Stokes  equations, 
explicit  or  implicit,  conservative  or  non-conservative,  with  chemistry  effects,  on  moving  or  stationary 
grids.  Numerous  methods  and  schemes  have  either  been  developed  or  applied  in  order  to  arrive  at  the 
solution  of  the  full  Navier-Stokes  equations.  These  methods  include,  but  are  in  no  way  limited  to  the 
handling  of:  turbulence,  chemistry,  Euler  fluxes,  time  integration,  convergence  acceleration,  matrix 
inversions  (implicit  solution),  and  zonal  interfaces. 

In  recent  years  the  NPARC  alliance  has  been  involved  with  the  development  of  the  WIND  code, 
a  finite  volume,  structured,  multi-zone,  compressible  flow  solver  with  flexible  chemistry  and  turbulence2 
The  code  incorporates  many  of  the  schemes  developed  over  the  years  which  have  demonstrated  varying 
degrees  of  accuracy  and  efficiency.  This  paper  presents  results  obtained  in  an  extensive  validation  effort 
using  many  combinations  of  the  different  schemes/methods  available  within  the  code  for  incompressible 
turbulent  flat  plate  flow.  The  numerical  results  are  compared  with  experimental  data  obtained  by 
Wieghardt3  for  friction  coefficient  and  boundary  layer  velocity  profiles.  The  results  show  that  most 
schemes  work  effectively,  and  that  efficiency  (in  terms  of  CPU  time)  varies  based  on  the  combination  of 
the  different  methods  applied  in  arriving  at  the  converged  solution. 

Test  Matrix 

Eight  key  factors  were  the  primary  focus  of  this  investigation  (WIND  defaults,  which  will  be 
referred  to  as  the  “baseline”  case,  are  depicted  using  [brackets]): 


1) .  y+  values  N/A 

2) .  Wall  functions  [OFF] 

3) .  Turbulence  models  [SST] 

4) .  Explicit  Euler  flux  schemes  [Roe  Physical] 

5) .  Euler  flux  orders  [2nd] 

6) .  Implicit  operators  [ADI] 

7) .  Time  integration  schemes  [Euler] 

8) .  Dimension  (2D  vs.  3D)  N/A 


The  strategy  of  this  investigation  will  be  categorized  as  follows: 


Case  1)  WIND  Version  2.55 

Starting  with  the  same  grid  used  in  previous  NPARC  validation  efforts  and  holding  the  number 
of  grid  points  constant,  y+  values  of  0.5,  1.0,  2.0,  4.0,  8.0,  16.0,  32.0,  64.0,  128.0,  256.0,  and  512.0  were 
investigated  for  the  “Baseline”  case:  wall  functions  [OFF],  turbulence  model  [SST],  Euler  flux  scheme 
[Roe  Physical],  order  [2nd],  implicit  operator  [ADI],  time  integration  scheme  [Euler],  and  dimension 
[2D]. 

Case  2)  WIND  Version  2.55 

Starting  with  the  same  grid  used  in  previous  NPARC  validation  efforts  and  holding  the  number 
of  grid  points  constant,  y+  values  of  0.5,  1.0,  2.0,  4.0,  8.0,  16.0,  32.0,  64.0,  128.0,  256.0,  and  512.0  were 
investigated  for  the  “Baseline”  case:  wall  functions  [ON],  turbulence  model  [SST],  Euler  flux  scheme 
[Roe  Physical],  order  [2nd],  implicit  operator  [ADI],  time  integration  scheme  [Euler],  and  dimension 
[2D]. 

Case  3)  WIND  Version  2.55 

Using  the  y+  =  1.0  grid  with  wall  functions  [OFF]  and  holding  the  explicit  Euler  flux  scheme 
[Roe  Physical],  order  [2nd],  implicit  operator  [ADI],  time  integration  scheme  [Euler],  and  dimension  [2D] 
fixed,  the  different  TURBULENCE  MODELS  were  studied,  namely: 

a) .  SST 

b) .  Spalart-Allmaras4 

c) .  Baldwin-Barth5 

d) .  PDT6 

e) .  Baldwin-Thomas7 

**Note:  Chien  k-s  was  previously  validated 

Case  4)  WIND  Version  2.43 

Using  the  y+  =  1.0  grid  with  wall  functions  [OFF]  and  holding  the  turbulence  model  [SST],  order 
[2nd],  implicit  operator  [ADI],  time  integration  scheme  [Euler],  and  dimension  [2D]  fixed,  the  different 
explicit  EULER  FLUX  SCHEMES  were  studied,  namely: 

a) .  Roe8 

b) .  van  Leer9 

c) .  Coakley10 

d) .  Roe  Physical" 

e) .  Central 

Case  5)  WIND  Version  2.55 
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iid  table 


function  iid_table  =  iid_table(n,  m,  vare,  number) 

%  function  iid_table  =  iid_table(n,  m,  vare,  number) 

"6 

%  Creates  number  test  matrices  which  are  n  X  m  with  mean  0,  variance 
%  var,  and  no  x  or  y  trend.  It  returns  a  number  X  4  matrix  where  the 
%  first  column  is  the  minimum  of  the  st  dev  statistic  for  the  test 
%  matrix  associated  with  the  row  number,  the  second  column  is  the 
%  mean  of  the  st  dev  statistic  for  the  test  matrix  associated  with 
%  the  row  number,  the  third  column  is  the  maximum  of  the  st  dev 
%  statistic  for  the  test  matrix  associated  with  the  row  number,  and 
%  the  fourth  column  is  the  standard  deviation  of  the  the  st  dev 
%  statistics  for  the  test  matrix  associated  with  the  row  number. 

for  i  =  1: number 

A  =  test_data (n,m,  0,  vare,  0,  0) ; 

X  =  sdev_stat_gen ( A) ; 

minn  =  min (X) ; 

meen  =  mean (X) ; 

maxx  =  max(X); 

vari  =  var (X) ; 

variance  =  var(vari); 

iid_table (i, 1)  =  min(minn); 

iid_table (i, 2)  =  mean(meen); 

iid_table (i, 3)  =  max (maxx); 

iid_table (i, 4 )  =  (variance) A . 5; 

end 
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Using  the  y+  =  1.0  grid  with  wall  functions  [OFF]  and  holding  the  turbulence  model  [SST], 
explicit  Euler  flux  scheme  [Roe],  implicit  operator  [ADI],  time  integration  scheme  [Euler],  and  dimension 
[2D]  fixed,  the  different  EULER  FLUX  ORDERS  were  studied,  namely; 


a) .  First-order 

b) .  Second-order 

c) .  Third-order 

d) .  Central 

e) .  Mixed  fully  upwind  and  central 

f) .  Third-order  fully  upwind 

g) .  Fourth-order  upwind-biased 

h) .  Fourth-order  central 

i) .  Fifth-order  upwind  biased 


Case  6)  WIND  Version  2.55 

Using  the  y+  =  1.0  grid  with  wall  functions  [OFF]  and  holding  the  turbulence  model  [SST], 
explicit  Euler  flux  scheme  [Roe  Physical],  order  [2nd],  implicit  operator  [ADI]  and  dimension  [2D]  fixed, 
the  different  TIME  INTEGRATION  SCHEMES  were  studied,  namely; 


a) .  ADI 

b) .  Gauss-Seidel  line  relaxation 

c) .  Jacobi  point  iterative  method 

Case  7)  WIND  Version  2.55 

Using  the  y+  =  1.0  grid  with  wall  functions  [OFF]  and  holding  the  turbulence  model  [SST], 
explicit  Euler  flux  scheme  [Roe  Physical],  order  [2nd],  implicit  operator  [ADI]  and  dimension  [2D]  fixed, 
the  different  TIME  INTEGRATION  SCHEMES  were  studied,  namely; 


a) .  Euler 

b) .  3  Stage  Runge-Kutta 

c) .  4  Stage  Runge-Kutta 

Case  8)  WIND  Version  2. 74 

For  the  “baseline”  case  a  3D  simulation  was  conducted.  The  3D  grid  consisted  of  five  planes 
spaced  Az  =  0.5”  apart,  where  each  plane  consisted  of  the  same  x  and  y  data  points  used  in  Case  1)  -  Case 
5).  The  3D  INDEX  SEQUENCING  was  investigated  in  order  to  determine  its’  effect  on  overall 
efficiency  and  performance  (CPU  time). 


a) .  ijk  sequencing 

b) .  jki  sequencing 

c) .  kij  sequencing 
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The  experimental  data  of  Wieghardt3,  namely  skin  friction  coefficient  and  boundary  layer 
velocity  profiles,  was  used  to  assess  the  numerical  results  of  the  validation  cases  presented  herein.  The 
flow  conditions  reported  by  Wieghardt,  and  used  in  the  numerical  simulations,  are  presented  in  Table  1. 


Table  8-1.  Flow  conditions 


Total  Pressure 

Static 

Total 

Static 

Mach  Number 

Reynolds 

(psia) 

Pressure 

Temperature 

Temperature 

Number 

(psia) 

(R) 

(R) 

(/ft) 

15.1 

14.7 

534.0 

530.0 

0.2 

1.38xl06 

Computational  Grids 

The  initial  2D  computational  grid  selected  for  this  investigation  was  the  same  as  that  used  in  the 
NPARC  validation  efforts,  and  is  depicted  in  Figure  8-1.  The  grid  consists  of  126  grid  points  in  the 
streamwise  direction  and  76  grid  points  in  the  transverse  direction.  The  grid  shown  in  Figure  8-1  is 
packed  near  the  wall,  with  approximately  5 1  points  placed  within  the  boundary  layer  while  maintaining  y+ 
=  1.0. 

A  grid  refinement  study  was  also  conducted  (Case  1).  Maintaining  the  same  number  of  grid 
points  in  the  streamwise  and  transverse  directions,  y+  values  of  0.5  -  512.0  (successively  doubling  y+,  see 
Case  1  description)  were  investigated. 

3D  simulations  were  also  conducted  for  y+  =  1.0.  The  3D  grid  consisted  of  five  planes  spaced  Az 
=  0.5”  apart,  where  each  plane  consisted  of  the  same  x  and  y  data  points  used  in  Case  1)  -  Case  5). 

Boundary  Conditions 

The  *.cgd  files,  which  contain  the  grid  data  and  boundary  condition  information,  were  created 
using  GRIDGEN  and  Gman.  Along  the  i  =  1  boundary,  the  “FROZEN”  boundary  condition  was  invoked 
and  the  Mach  number  (0.2),  total  pressure  (15.1 157  psia),  total  temperature  (523.152  R),  angle  of  attach 
(0),  and  side  slip  angle  (0),  were  specified.  Along  the  i  =  imax  boundary,  the  “CONFINED  OUTFLOW” 
boundary  condition  was  implemented.  Along  the  j  =  1  boundary,  the  “INVISCID  WALL”  condition  was 
invoked  for  i  =  1  -  25,  whereas  the  “VISCOUS  WALL”  condition  was  specified  for  i  =  26  -  126.  Finally, 
along  the  j  =  jmax  boundary,  the  “INFLOW/OUTFLOW”  condition  was  used. 

For  the  3D  case  ijk  sequencing,  the  above  is  still  true  where  the  i  and  j  lines  are  now  i  and  j 
constant  planes.  For  the  k  =  1  plane  and  k  =  kmax  plane,  the  reflection  boundary  condition  was  specified. 

Convergence 

For  each  case  investigated,  a  CFD  number  of  10  was  initially  chosen.  The  solution  was  advanced 
and  considered  converged  when  the  residual  began  to  plateau  (see  Figure  8-2).  The  CFL  condition  was 


1-6 


than  incremented  by  10,  the  flowfield  reinitialized  to  freestream  conditions,  and  the  solution  advanced.  If 
the  solution  proved  stable,  the  CFL  number  was  then  incremented  further  and  the  process  repeated  until 
the  maximum  CFL  number  was  determined.  Thus,  the  maximum  allowable  CFL  condition  for  each  case 
was  determined,  and  convergence  was  guaranteed.  It  was  determined  that  only  the  RHS  operator  had  an 
effect  on  the  maximum  CFL  number  that  could  be  used.  It  should  be  noted  that  only  the  Roe  and  Roe 
physical  schemes  allowed  CFL  numbers  greater  than  10,  in  fact  one  order  magnitude  higher  than  the 
other  schemes  investigated. 

Results 

Figure  8-3  shows  the  predicted  friction  coefficient  and  velocity  profiles  at  5  “i”  locations  obtained 
using  both  the  NPARC  and  WIND  codes.  In  this  particular  case,  the  explicit  flux  operator  [CENTRAL], 
implicit  operator  [ADI],  Euler  flux  order  [2nd],  and  grid  [i  =  256,  j  =  76,  and  y+  =  1.0]  are  the  same,  and 
the  turbulence  models  are  similar  [SST  for  WIND  (which  is  a  mixed  k-w  formulation  in  the  boundary 
layer  and  k-e  formulation  in  the  freestream)  and  Chein  k-e  model  for  NPARC].  In  terms  of  the  friction 
coefficient,  the  figure  shows  that  the  results  obtained  with  WIND  are  in  much  closer  agreement  with  the 
experimental  data  than  was  obtained  using  the  NPARC  code.  However,  at  the  i  =  32  location,  which  is 
near  the  plate  leading  edge,  the  opposite  is  true.  Note  that  as  the  velocity  profile  develops,  the  numerical 
results  obtained  using  the  two  codes  are  nearly  indistinguishable,  and  both  sets  of  results  are  in  very  good 
agreement  with  the  experimental  data. 

Case  1)  Figure  8^f  shows  the  numerical  results  obtained  using  the  WIND  code  in  terms  of  the 
friction  coefficient  and  velocity  profiles  at  5  “i”  locations  for  the  different  y+  grids.  It  can  be  seen  that  y+ 
values  up  to  about  4.0  yield  good  results  when  compared  to  the  experimental  data.  As  the  value  of  y+ 
increases,  the  first  grid  point  near  the  wall  moves  out  of  the  viscous  sublayer  and  the  assumptions  the 
various  transport  turbulence  models  are  based  on  become  invalid.  Figure  8-4  also  shows  the  variation  of 
y+  with  Reynold’s  number.  Notice  that  y+  does  vary  somewhat  from  the  plate  leading  edge  to  the  trailing 
edge. 

Case  2)  Figure  8-5  compares  the  effect  of  y+  with  and  without  wall  functions.  As  can  be  seen  in 
the  figure,  the  wall  function  allows  for  much  larger  values  of  y+  while  still  yielding  fairly  accurate  results. 
Note  that  near  the  plate  leading  edge,  a  kink  is  evident.  For  y+  =  256.0  it  occurs  at  Re  =  0.8  x  106 ,  at  y+  = 
64.0  it  occurs  at  Re  =  '1.2  x  106  ,  and  at  y+  =  32.0  it  occurs  at  Re  =  2.8  x  106.  This  can  most  likely  be 
attributed  to  the  fact  that  y+  is  not  constant  along  the  plate  surface,  as  was  depicted  in  Figure  8-5.  As  y+ 
increases  with  Reynolds  number,  points  are  pulled  out  of  the  viscous  sublayer  and  the  wall  functions  at 
which  point  the  wall  functions  are  switched  on. 

Case  3)  Figure  8-6  displays  the  results  obtained  using  different  turbulence  models  while 
maintaining  everything  else  the  same.  Excellent  agreement  among  most  of  the  different  turbulence 
models  and  Wieghardt’s  data  can  be  observed  in  terms  of  the  friction  coefficient.  There  is  also  good 
agreement  in  terms  of  the  velocity  profile  except,  perhaps,  for  i  =  32,  near  the  plate  leading  edge.  There 
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does  appear  to  be  a  problem  with  the  Baldwin-Barth  model.  It  severely  underpredicts  the  friction 
coefficient  and  seems  to  lack  sufficient  dissipation  near  the  outer  edge  of  the  boundary  layer.  This  is  an 
indication  of  a  possible  bug  in  the  code. 

Case  4)  Figure  8-7  presents  the  results  obtained  using  the  different  RHS  operators.  Again, 
excellent  agreement  between  the  numerical  results  and  the  experimental  data  is  observed  in  terms  of  the 
friction  coefficient  (note  that  the  van  Leer  data  is  not  seen  in  this  plot).  It  seems  that  the  only  discrepancy 
exists  near  the  plate  leading  edge,  as  was  observed  in  Case  5.  While  the  central  scheme  is  too  dissipative, 
the  Roe  scheme  is  lacking  in  sufficient  dissipation.  Note  that  the  van  Leer  schemes  produces  poor  results 
at  each  i  location,  indicating  a  possible  bug  in  the  code. 

Case  5)  Figure  8-8  presents  the  results  obtained  using  the  different  RHS  orders.  Note  that  it  was 
determined  that  the  third  order  fully  upwind  order  was  not  functioning  correctly.  The  code  terminated  the 
first  iteration  for  no  apparent  reason.  Other  than  3U,  all  other  orders  produced  excellent  results  as  can  be 
seen  in  the  figure. 

Case  6)  In  this  case,  not  only  was  the  accuracy  of  interest,  but  the  time  required  for  convergence 
as  well.  Figure  8-9  presents  the  results  obtained  using  the  different  implicit  operators.  As  can  be  seen  in 
Figure  8-9,  there  is  little,  if  any,  discernible  difference  between  the  different  numerical  results.  However, 
as  can  be  seen  in  Figures  8-10,  the  time  required  for  such  solutions  varies  greatly.  While  Figure  8-10 
assures  that  the  convergence  rate  is  identical  for  the  3  different  implicit  operators,  the  time  required  for 
convergence  is  much  lower  for  the  ADI  scheme  than  the  2  iterative  methods.  Note,  furthermore,  that  the 
approximate  factorization  error  of  the  ADI  scheme,  does  not  have  a  noticeable  effect  on  any  of  the  results 
obtained  for  this  particular  case. 

Case  7)  As  with  Case  6,  this  case  was  also  concerned  with  the  efficiency  of  the  different 
numerical  methods  as  well  with  the  accuracy  of  the  results.  Figure  8-11  shows  that  no  discernible 
difference  exists  in  terms  of  the  friction  coefficient  and  velocity  profiles.  However,  it  was  determined  that 
the  Euler  scheme  proved  to  be  more  than  50%  faster  than  the  3  stage  Runge  Kutta  scheme  and  nearly  70% 
faster  than  the  4  stage  Runge  Kutta  scheme.  For  a  steady  state  problem,  it  is  quite  obvious  that  the  order 
accuracy  is  not  nearly  as  important  as  the  time  require  for  convergence. 

Case  8)  As  a  final  test  of  the  WIND  code,  some  3D  simulations  were  investigated  which  yielded 
interesting  results.  First,  the  sequencing  of  the  indicies  was  varied  to  determine  which  resulted  in  the 
fastest  convergence.  Figure  8-12  shows  that,  in  terms  of  the  friction  coefficient,  there  exists  no  difference 
between  the  different  planes  and  different  sequencing.  Note,  however,  that  the  3D  friction  coefficient  is 
higher  than  the  2D  friction  coefficient,  keeping  all  numerical  schemes/methods  the  same.  Also,  note  the 
“kink”  that  exists  near  Re  =  5.0  x  106.  Figure  8-13  presents  the  CPU  requirements  for  the  different 
sequencing.  It  can  be  seen  that  the  jki  sequencing  is  approximately  3%  faster  than  the  ijk  sequencing  and 
nearly  6%  faster  than  the  kij  sequencing.  Figure  8-14  presents  the  convergence  rates  for  the  different 
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sequencing,  allowing  one  to  conclude  that  the  same  number  of  iterations  are  required  for  convergence  and 
that  the  results  presented  in  Figure  8-13  is  in  fact  an  accurate  assessment  of  efficiency. 

It  was  hypothesized  that  the  “kink”  seen  in  Figure  8-12  could  possibly  be  attributed  to  the 
turbulence  model  (SST).  As  such,  a  case  was  also  investigated  using  the  PDT  and  Spalart-Allmaras 
turbulence  models,  and  the  3D  results  compared  with  the  2D  results.  First,  in  looking  at  Figure  8-15,  note 
that  convergence  appears  to  occur  much  faster  using  the  algebraic  model  than  the  2  transport  models,  and 
that  the  1  equation  model  appears  to  converge  slightly  faster  than  the  2  equation  model,  as  one  would 
expect.  However,  looking  at  the  results  presented  in  Figure  8-16,  it  can  be  seen  that  the  Spalart-Allmaras 
model  yields  the  best  agreement  when  comparing  2D  and  3D  results.  Also,  the  “kink”  is  no  longer 
apparent.  Finally,  it  seems  that  there  is  a  problem  with  the  PDT  model  in  that  the  friction  coefficient 
appears  almost  laminar.  There  is,  at  present,  no  explanation  why  the  2D  and  3D  simulations  differ  when 
changing  the  different  turbulence  models. 

Conclusions 

A  validation  effort  using  the  WIND  code  to  solve  for  incompressible  turbulent  flow  past  a  flat 
plate  was  conducted.  The  results  show  that  most  numerical  methods  and  schemes  work  effectively  with 
the  exception  of:  the  Baldwin-Barth  turbulence  model,  the  3rd  Order  fully  upwind  order  for  the  RHS 
operator,  and  the  van  Leer  RHS  operator.  It  was  observed  that  efficiency  (defined  as  the  number  of 
iterations  per  CPU  second)  more  than  accuracy  varies  depending  on  the  combination  of  the  different 
schemes/methods  for  this  simple  test  case.  It  would  be  interesting  to  conduct  this  same  investigation, 
using  the  same  test  matrix,  for  a  slightly  more  complicated  test  case  where  shocks,  shear  layers  and 
separation  are  present.  Only  than  can  one  conclude  the  best  combination  of  schemes/methods  to  be  used 
for  a  particular  case. 
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ABSTRACT 

Proper  functioning  of  the  mammalian  nervous  system  requires  myelination  of  neuronal  axons. 
Myelination  of  the  newborn  mouse  brain  begins  shortly  after  birth  and  is  complete  by  about  20  days 
of  age.  Abnormalities  of  myelin  basic  protein  (MBP)  production  have  a  direct  impact  on  myelination 
as  has  been  demonstrated  in  the  shiverer  mouse.  Mice  with  less  than  25%  of  normal  MBP  levels  have 
aberrant  myelination  and  brain  development  and  demonstrate  a  characteristic  tremor  at  12  days  of 
age.  Growth  hormone  (GH)  and  insulin-like  growth  factor  I  (IGF-1)  have  been  shown  in  vitro  and 
in  vivo  to  affect  myelination.  This  project  explores  the  effect  of  GH  deficiency  on  myelin  basic 
protein  expression  to  determine  whether  GH  deficiency  exacerbates  MBP  haploinsufficiency.  These 
studies  grew  out  of  studies  undertaken  on  children  with  1 8q-  syndrome  who  have  only  a  single  copy 
of  the  MBP  gene,  are  hypomyelinated,  and  are  also  shown  to  suffer  growth  hormone  deficiency  or 
insufficiency.  A  hybrid  mouse  model  that  mimics  these  deficiencies  of  1 8q-  patients  was  developed 
to  explore  the  relationship  between  GH  and  MBP  in  vivo. 
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INTRODUCTION 

Expression  of  the  myelin  basic  protein  (MBP)  gene  is  critical  for  normal  myelination  of  the  central 
nervous  system  (CNS).  This  has  been  most  clearly  demonstrated  by  the  neuropathology  of  the 
mutant  shiverer  ( shi/shi )  mouse  in  which  a  large  portion  of  the  MBP  gene  is  deleted  and  MBP  protein 
is  essentially  undetectable  (1,2).  Consequently,  the  CNS  has  virtually  no  compact  myelin.  These 
mutants  exhibit  a  characteristic  tremor  (shiver)  appearing  at  12  days  of  age,  coincident  with  the 
normal  age  at  which  myelination  occurs  in  rodents,  and  die  prematurely  usually  during  status 
epilepticus.  That  MBP  deficiency  is  responsible  for  the  phenotype  is  supported  by  mutant  rescue  with 
a  MBP  transgene  (3). 

Cross-breeding  experiments  utilizing  mice  transgenic  for  the  MBP  gene  have  established  the  level  of 
MBP  expression  sufiicient  for  normal  myelinogenesis  (4).  Heterozygous  shiverer  mice  produce  about 
50%  of  the  normal  level  of  MBP  mRNA  and  a  proportionate  amount  of  MBP  protein.  These  animals 
exhibit  a  normal  behavioral  and  morphological  phenotype;  however,  they  have  minor  biochemical 
changes  in  myelin  composition.  Mice  which  express  less  than  25%  of  MBP  display  the  shivering  of 
the  dysmyelinating  phenotype.  These  data  indicate  that  the  level  of  MBP  expression  influences  the 
assembly  of  myelin  by  oligodendrocytes  which  suggests  that  factors  which  influence  MBP  expression 
also  influence  myelin  synthesis  and  assembly.  Additionally,  factors  that  regulate  oligodendroglial 
proliferation  and  differentiation  may  affect  myelinogenesis.  In  this  regard,  there  is  much  evidence, 
both  in  vivo  and  in  vitro,  that  growth  hormone  (GH)  and  insulin-like  growth  factor  1  (IGF-1)  enhance 
myelination  (5-10).  GH-deficient  mice,  such  as  the  little  (lit/lit)  mutant  which  has  an  abnormality  of 
the  GH-releasing  hormone  receptor  (Ghrhr)  (11),  have  been  reported  to  be  hypomyelinated  and 
correctable  by  early  postnatal  GH  treatment  (12). 

Based  on  these  observations,  GH  may  be  a  potentiator  of  MBP  production.  To  test  this,  we  have 
examined  the  impact  of  naturally  occurring  GH  deficiency  ( liVlit )  on  MBP  and  myelination  in  vivo, 
investigated  the  effects  of  MBP  deficiency  (shi/shi)  on  GH  and  IGF-1  production,  and  explored  the 
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relationship  between  GH  and  MBP  in  a  hybrid  mouse  (lit/lit,  -shi)  model.  We  hypothesized  that 
if  GH  affects  MBP  production,  then  MBP  levels  that  are  already  reduced  in  the  +/shi  mouse  may  be 
further  reduced  in  the  GH-deficient  hybrid,  and  therefore,  the  affected  mice  may  display  the 
dysmyelinating  phenotype.  In  this  case,  early  treatment  with  GH  might  delay  or  prevent  the 
appearance  of  the  tremor.  Alterations  in  developmental  progression  of  myelination  and  its  impact  on 
function  are  being  investigated. 

These  studies  in  mice  grew  out  of  studies  undertaken  on  children  with  lBq-  syndrome  (1 1,12,13). 
The  characteristic  features  of  this  syndrome  are  short  stature,  mental  retardation  and  deafness. 
Affected  children  have  only  a  single  copy  of  the  MBP  gene  and  are  hypomyelinated  as  demonstrated 
by  quantitative  magnetic  resonance  relaxometry  (16).  Half  of  affected  children  are  GH-deficient  and 
the  other  children  have  evidence  of  dysregulation  of  GH  production.  Accumulating  human  data 
suggest  that  GH  treatment  may  improve  both  myelination  and  cognitive  function  (unpublished)  raising 
the  possibility  that  a  relationship  exists  between  MBP  and  GH.  These  mouse  studies  permit  the 
opportunity  to  more  thoroughly  explore  the  relationship  between  GH  and  MBP  within  the  context 
of  the  developing  brain  and  may  provide  insight  into  novel  treatment  approaches  to  children  with 
hypomyelination. 

METHODOLOGY 

Mouse  colony:  Little  and  shiverer  mutant  mice  were  supplied  by  Jackson  Laboratory.  The 
little  mutant  is  in  the  C57B1/6J  background,  and  the  shiverer  mutant  is  in  the  C3HeB/FeJ  background. 
Wild-type  little  males  were  crossed  with  lit/lit  females,  and  the  heterozygous  offspring  (FI)  were 
interbred  to  create  +/+,  lit/+,  and  lit/lit  genotypes.  C3H  shiverer  mice  were  also  interbred  to 
produce  C3H  mice  for  study.  The  hybrid  mouse  (lit/lit; shi/ +)  was  produced  by  first  crossing  C57B1 
lit/ lit  females  with  C3H  shi/ shi  males  to  produce  a  double  heterozygote.  Male  FI  mice  were 
backcrossed  to  C57B1  lit/lit  females.  Mice  were  housed  in  the  animal  facility  at  the  Clinical 
Investigations  Directorate  of  Wilford  Hall,  Lackland  Air  Force  Base,  fed  standard  laboratory  chow 
ad  libitum  and  kept  on  a  12  h  dark/light  schedule. 

Mutant  genotyping:  A  method  was  developed  to  genotype  the  little  locus  (Ghrhr)  by  PCR  and 
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subsequent  restriction  analysis  (17).  The  shiverer  locus  (MBP)  was  determined  by  Southern  analysis 
as  previously  published  (18).  A  581 -bp  probe  that  spanned  the  deletion  breakpoint  was  designed  for 

this  purpose. 

Evaluation  of  brain  MBP  and  PLP  mRNA  and  MBP  protein :  Brains  were  harvested  immediately 
after  death  induced  by  carbon  dioxide  inhalation.  The  left  hemisphere  was  frozen  in  liquid  nitrogen 
and  stored  at  -20°C  for  protein  analysis  and  the  right  hemisphere  was  homogenized  in  4M 
guanidinium  thiocyanate  buffer  for  RNA  isolation.  Total  RNA  was  extracted  using  a  protocol 
adapted  from  Chomczynski  and  Sacchi  (19).  Five  pg  of  total  RNA  were  electrophoresed  through 
a  1.2%  denaturing  agarose  gel,  transferred  onto  a  GeneScreen  membrane  (DuPont  NEN)  and  UV 
crosslinked.  The  membrane  was  hybridized  with  a  radiolabeled  cDNA  probe  for  the  14kD  form  of 
mouse  MBP  (gift  of  Dr.  Carol  Readhead,  Cedar  Sinai  School  of  Medicine)  according  to  Church  and 
Gilbert  (20).  The  membrane  was  exposed  for  24  hr  to  a  phosphorimager  screen  and  quantitation  was 
done  by  ImageQuant  version  3.2  (Molecular  Dynamics,  CA).  The  membrane  was  subsequently 
stripped  in  50mM  sodium  phosphate  in  55%  formamide  at  65°C  for  60  min.  and  rehybridized  with 
a  radiolabeled  cDNA  probe  for  mouse  cyclophilin  (Ambion,  Inc.,  Austin  TX).  Following  quantitation 
and  stripping,  the  membrane  was  hybridized  with  a  radiolabeled  cDNA  probe,  BAS  101 3,  for  mouse 
PLP  (gift  of  Dr.  Anthony  Campagnoni,  UCLA).  Since  it  has  been  reported  (21)  that  cyclophilin 
mRNA  levels  slightly  decrease  over  the  first  40  days  of  postnatal  life,  the  membranes  were  also 
rehybridized  with  a  cDNA  probe  for  a-tubulin  (21).  The  message  levels  for  MBP  and  PLP  were 
normalized  to  the  cyclophilin  mRNA  abundance  or  to  the  a-tubulin  mRNA  abundance.  Since  the 
relative  values  were  the  same,  only  the  message  levels  normalized  to  cyclophilin  abundance  are 
reported. 

MBP  protein  analysis  was  done  using  a  radioimmunoassay  kit  purchased  from  Diagnostic 
Systems  Laboratories  (Webster  TX).  The  frozen  left  hemisphere  of  each  brain  was  thawed  on  ice 
and  immediately  homogenized  in  5  ml  cold  sterile  H20.  The  homogenate  was  centrifuged  at  3000g 
for  20  min  at  4°C.  The  supernatant  was  diluted  in  0.2M  Tris-acetate  buffer,  pH  7.2,  containing  1% 
Triton  X-100  and  0.1%  aprotinin  (22),  and  an  aliquot  was  used  to  measure  MBP  concentration. 
MBP  concentration  was  normalized  to  total  protein  concentration  as  measured  by  the  Bio-Rad  DC 
colorimetric  assay  (Richmond  CA). 
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IGF- 1,  GH  and  TSH  quantitation :  Blood  was  collected  from  mice  by  cardiac  puncture. 
Serum  was  separated  from  the  cellular  components  and  frozen.  Fifty  p.1  of  serum  was  used  per  assay; 
assays  were  performed  in  duplicate.  Each  sample  was  acid-extracted  using  C-l  8  cartridges  (Sep-Pak, 
Waters  Associates)  and  IGF-1  assayed  by  radioimmunoassay.  Serum  assay  of  IGF-1  was  performed 
by  Dr.  Ross  Clark  at  Genentech,  Inc.,  San  Francisco,  CA.  Anterior  pituitaries  were  harvested  for 
analysis  of  intrapituitary  GH  and  TSH,  weighed,  homogenized  in  500  (il  0.05M  PBS  with  1%  BSA, 
and  stored  at  -20°C.  At  the  time  of  assay,  the  homogenate  was  centrifuged  at  8,000g  for  20  min  at 
4°C.  The  supernatant  was  used  for  analysis  with  enzyme  immunosorbent  assay  (EIA)  kits  purchased 
from  Amersham  (Arlington  Heights,  IL)  for  either  rat  GH  or  TSH. 

Histological  analysis  of  CNS  myelin:  60-day  old  mice  were  euthanized  by  isoflurane 
inhalation,  and  brain  and  spinal  cord  were  removed  and  fixed  in  10%  buffered  formalin  for  1-3  days. 
The  brains  were  then  bisected  and  placed  into  histology  cassettes.  The  spinal  cord  tissue  was 
bisected  to  show  a  cross  section  and  sectioned  lengthwise  depending  on  the  size  of  the  tissue.  The 
tissues  were  processed  in  a  Shandon  Lipshaw  processor  according  to  standard  techniques  (Carson, 
1 997).  They  were  than  embedded  in  paraffin  and  sectioned.  Semithin  sections  of  8(im  of  the  brains 
were  cut  sagittally,  and  the  spinal  cord  was  sliced  at  4pm.  The  sections  were  then  deparaffinized  and 
hydrated  to  95%  ethanol  and  placed  in  0.1%  luxol  fast  blue  solution  overnight  at  55°C.  The  slides 
were  rinsed  in  95%  ethanol  followed  by  distilled  water  and  then  dipped  in  lithium  carbonate  and 
differentiated  in  70%  ethanol.  They  were  counterstained  with  cresyl  echt  violet.  Size  measurements 
of  white  matter  regions  were  made  with  the  assistance  of  ImagePro  (Media  Cybernetics,  Silver 
Springs,  MD). 

Electron  microscopy  of  optic  nerves  and  brain:  60-day  mice  and  controls  (n=5  for 
each)  were  deeply  anesthetized  with  sodium  pentobarbital  (100  mg/kg),  and  transcardically  perfused 
with  1%  paraformaldehyde  and  2.5%  glutaraldehyde  in  0. 1M  phosphate  buffer,  ph  7.4,  at  25°C.  The 
brain  and  the  optic  nerve  connected  to  the  eye  were  removed  and  fixed  overnight  in  the  same  fixative 
at  4°C.  The  brain  was  sagittally  sliced  in  ~0.5  mm  of  thickness  in  the  region  containing  the  anterior 
commissure  and  the  corpus  callosum.  The  tissues  were  then  washed  in  0.1M  phosphate  buffer, 
postfixed  in  1%  OsCL  in  the  same  buffer,  rinsed  with  H2O,  dehydrated  through  a  graded  series  of 
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ethanol  and  propylene  oxide,  and  directionally  embedded  in  epoxy  resin  (Polybed  812,  Polysciences, 
Inc,  Warrington,  PA).  Semithin  sections  of-l^im  of  the  optic  nerves  were  cut  perpendicular  to  the 
longitudinal  axis  of  the  nerves  at  a  point  1  mm  behind  the  retina.  Sections  were  stained  with  7% 
aqueous  uranyl  acetate  and  counterstained  with  Reynolds’  lead  citrate  (Hayat,  1970)  using  a 
microwave  staining  procedure  (25).  Semithin  sections  of  the  brain  tissue  containing  the  anterior 
commissure  were  cut  sagittally  and  processed  in  the  same  manner.  The  sections  were  examined  with 
a  Philips  301  transmission  electron  microscope  (Philips  Electronic  Instruments  Co,  Mahwah,  NJ, 
USA). 

Six  to  eight  photographs  were  randomly  taken  from  each  anterior  and  posterior  part  of  the 
anterior  commissure  and  from  each  corpus  callosum  at  10,000X.  Each  negative  was  printed  to  give 
a  final  30,000X  magnification.  The  axon  size  and  thickness  of  myelin  sheaths  were  measured  with 
the  aid  of  the  Image-Pro  analysis  system.  Only  those  axons  that  were  perpendicular  to  their  long  axis 
■  and  in  which  both  the  inner  and  outer  surface  of  the  myelin  sheath  was  sharp  and  distinct  were  used 
in  this  analysis.  The  diameter  of  axons  was  calculated  from  the  circumference  of  the  inner  surface 
of  the  myelin  sheaths.  The  thickness  of  myelin  sheaths  was  measured  by  determining  the 
circumference  of  the  outer  surface  of  the  sheath  and  using  the  approximation  as  calculated  by 
ImagePro  which  subtracts  the  axon  circumference  and  corrects  for  the  shape  of  the  image. 

Visual  Evoked  Potentials:  Flash  visual  evoked  potentials  (FVEPs)  were  recorded  from  the 
mutant  and  wild-type  control  (n=l  1  for  each)  mice.  The  mice  were  anesthetized  with  ketamine  and 
xylazine  (75  mg/kg  and  lOmg/kg,  respectively  for  wild-type  and  little  mice,  and  60  mg/kg  and  8 
mg/kg,  respectively  for  shiverers)  under  normal  room  illumination.  They  were  placed  prone  on  a 
sling,  and  normal  body  temperature  was  maintained  by  use  of  a  surrounding  hot  water  bag.  The  mice 
were  placed  facing  the  flash  lamp  diffusing  face  plate  (Grass  PS  22  Photostimulator)  at  a  distance  of 
20  cm.  Corneas  were  kept  moist  with  saline  or  2.5%  hydroxypropyl  methylcellulose  (Goniosol, 
Johnson  &  Johnson).  Stainless  steel  needle  electrodes  were  placed  subdermally  in  the  occipital 
midline  (active)  and  the  ear  (reference).  The  active  electrode  was  shortened  to  6  mm  length  and  re¬ 
sharpened  and  then  inserted  posteriorly  so  that  its  tip  was  not  anterior  to  the  interaural  line  in  an 
attempt  to  reduce  interference  by  the  electroretinogram  (ERG).  A  subdermal  needle  electrode  in  the 
midline  near  the  tail  served  as  the  ground  electrode.  The  signals  were  amplified  by  10,000  and  filtered 
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from  l-1000Hz  (Grass  P5  amplifier).  The  amplified  signals  containing  both  electroencephalographic 
and  electrocardiographic  activity  were  monitored  on  the  screen  of  an  oscilloscope  (Tektronix  5110). 
The  amplified  signals  were  also  digitized  and  recorded  by  a  Nicolet  4094  digital  oscilloscope  set  for 
averaging  (0.396  s  trace,  digitization  rate  =  10,000/s).  During  recording  of  FVEPs,  ambient 
illuminance  was  dim  (2.9  lx)  as  measured  at  the  eye  with  a  light  meter  (International  Light  IL  1700 
radiometer  with  illuminance  detector).  The  surrounding  luminance  was  0.3  to  0.6  fL  (Minolta  1 
degree  Spotmeter).  The  mice  were  stimulated  binocularly  at  1  flash/s.  Flash  intensity  was  set  at  16, 
which  produced  88  be  s  (Gossen  Luna  Pro  light  meter)  illumination  at  the  eyes.  One  hundred  twenty 
traces  synchronized  with  flash  onset  were  averaged.  Four  FVEPs  were  recorded  from  each  mouse. 
Noise  controls  were  produced  in  the  same  way  but  with  the  flash  occluded  by  a  black  opaque  cloth. 
The  digitized  waveforms  were  saved  as  ASCII  files  using  a  commercial  software  analysis  package 
(Vu-Point).  The  average  of  the  four  FVEPs  was  computed  and  used  as  the  measured  datum  for  each 
’  mouse. 

RESULTS 

Myelin  characterization  of  little  mice: 

MBP  and  PLP  levels  in  brain  from  little  mice 

Although  there  is  extensive  data  on  the  effect  of  GH/IGF-1  on  myelination,  there  is  no  information 
on  the  specific  effect  of  GH/IGF-1  on  MBP  in  vivo.  The  period  of  maximal  MBP  expression  in  mice 
is  between  10  and  30  days  after  birth  with  the  peak  at  about  15  days  (26).  Thereafter,  expression 
declines  and  plateaus  at  a  much  lower  level.  Therefore,  MBP  expression  was  determined  at  18,  25 
and  60  days  postnatal  in  little  mice.  As  seen  in  figure  2,  MBP  mRNA  expression  in  little  mice  is 
lower  than  wild-type  mice,  but  is  only  significant  at  60-days.  However,  there  is  no  difference  in  MBP 
protein  levels  between  little  and  wild-type  mice  at  any  age  tested.  Therefore,  any  hypomyelination 
of  little  mice  cannot  be  attributed  to  insufficient  MBP  expression. 
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Figure  1.  Northern  analysis  of  MBP  and  PLP  mRNA 
in  little  mice. 


MBP  expression  in  littie  mice 
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Figure  2.  Quantitation  of  MBP  mRNA  and  protein  in 
little  mice.  MBP  mRNA  data  are  shown  as  a  ratio  of 
MBP  mRNA  to  cyclophilin  mRNA. 


Since  MBP  protein  levels  were  found  to  be  normal  in  little  mice,  expression  levels  of  PLP  were  also 
determined.  Figure  3  shows  that  PLP  mRNA  levels  in  little  mice  were  no  different  from  wild-type; 
therefore,  protein  levels  were  not  measured.  These  results  indicate  that  low  PLP  is  not  responsible 
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for  any  hypomyelination  in  the  little  mice. 


PLP  tnPHA 


*  p  =  .0677 

Figure  3.  Quantitation  of  PLP  mRNA  in  little  mice. 

Data  are  shown  as  a  ratio  of  PLP  mRNA  to 
cyclophilin  mRNA. 

Morphometric  analysis  of  white  matter  tracts  and  optic  nerve  in  little  mice 
Our  molecular  analyses  of  myelin  proteins  are  inconsistent  with  the  literature  reporting 
hypomyelination  in  the  little  mice;  therefore,  a  morphometric  analysis  of  the  major  CNS  myelinated 
regions  was  done  in  order  to  confirm  or  refute  these  reports.  The  thickness  of  myelin  sheaths,  axon 
sizes,  and  density  of  myelinated  axons  were  measured  in  the  corpus  callosum,  anterior  commissure, 
and  the  optic  nerve  in  little  and  wild-type  mice.  The  corpus  callosum  and  anterior  commissure  were 
selected  for  study  because  these  regions  have  previously  been  shown  to  be  abnormally  myelinated  in 
IGF-1  knockout  and  IGFBP-1  mice  (8,9).  As  shown  in  tables  1  and  2,  our  results  indicate  no 
difference  in  either  the  thickness  of  myelin  sheaths  or  the  density  of  myelinated  axons  in  the  corpus 
callosum.  Since  myelin  thickness  increases  with  increasing  axon  diameter,  the  myelin  sheaths  of 
comparable  axon  sizes  were  compared.  Additionally,  no  difference  in  average  axon  diameter  in  these 
white  matter  regions  was  seen.  The  results  for  the  anterior  commissure  are  shown  in  tables  1  and  3. 
The  density  of  the  myelinated  axons  is  similar  between  little  mice  and  wild  controls  in  both  the 
anterior  and  the  posterior  region.  Moreover,  the  thickness  of  the  myelin  sheaths  in  both  regions  is 
not  reduced  in  the  little  mice.  In  contrast,  our  data  indicates  a  statistically  significant  increase  in 
average  myelin  thickness  of  axons  less  than  1000  nm  in  diameter  within  the  posterior  region  of  the 
anterior  commissure  of  the  little  mouse  as  compared  to  the  wild-type  controls. 
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Table  1.  The  density  of  myelinated  axons  in  little  and  wild-type  mice 


Corpus  callosum 

Anterior 

Anterior 

commissure 

commissure 

(anterior) 

(posterior) 

Little 

38.2  ±9.3 

64.4  ±26.6 

38.3  ±11.7 

Wild 

41.9  ±  12.5 

69.8  ±  14.9 

46.4  ±20.3 

The  number  of  myelinated  axons  was  counted  on  each  electron  microphotograph.  Values  are  means  ±  SD. 


Table  2.  Myelin  sheath  thickness  in  corpus  callosum 


<  500  nm 

500-800  nm 

>  800  nm 

Little 

Avg.  myelin 
thickness 
(Median) 

94.9  ±  15.3 

(93.9) 

105.6  ±25.0 

(101.1) 

131.3  ±41.7 

(119.4) 

Avg.  axon 
diameter 

399.8  ±71.2 
F23.4%1 

637.6  ±88.0 
F45.4%1 

1093.1  ±307.1 
F31.2%1 

Avg.  myelin 
thickness 

92.7  ±  15.4 

108.8  ±72.6 

131.4  ±  32.0 

Wild 

(Median) 

(90.7) 

(102.0) 

(126.8) 

Avg.  axon 

414.2  ±71.6 

641.2  ±82.0 

1027.5  ±233.8 

diameter 

T24.0%1 

r49.0%l 

F27.0%] 

Values  are  means  ±  SD.  The  percentage  of  axons  analyzed  is  shown  in  brackets. 

Table  3.  Myelin  sheath  thickness  in  anterior  commissure 

Anterior 

region 

<  500  nm 

500-700  nm 

>  700  nm 

Avg.  myelin 
thickness 

95.2  ±27.9 

105.9  ±25.2 

128.6  ±34.6 

Little 

(Median) 

(90.3) 

(103.4) 

(125.8) 

Avg.  axon 

410.8  ±75.2 

597.7  ±57.7 

974.4  ±274.6 

diameter 

F20.4%] 

F36.1%] 

**  F43 . 5%1 

Avg.  myelin 
thickness 

93.7  ±  21.6 

101.6  ±  27.3 

124.2  ±29.4 

Wild 

(Median) 

(89.5) 

(98.1) 

(120.4) 

Avg.  axon 

396  2  ±69.5 

584.8  ±58.5 

884.4  ±  169.6 

diameter 

F40.2%1 

F36.6%] 

F23.2%1 
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Posterior  region 


Avg.  myelin 
thickness 
Median) 
Avg.  axon 
diameter 


Avg.myelin 
thickness 
(Median) 
Avg.  axon 
diameter 


<  700  nm  I  700-1000  nm 


107.3  ±20.5  *  124. 2  ±25.1* 


(106.4) 
552.3  ±  106.8 
[25.2%] 


85.2  ±  19.7 
(84.0) 

543.2  ±  124.9 
37.7%] 


(1189) 
857.1  ±79.7 
[31.1%] 


107.2  ±23.9 

(103.2) 
845.0  ±86.3 
[39.4%] 


>  1000  nm 


144.0  ±31.3 
(144,7) 

1318.9  ±  273.3 
[43.7%] 


134.2  ±37.0 
(126.2) 

1228.8  ±210.7 

[22.8%] 


Values  are  means  ±  SD.  The  percentage  of  axons  analyzed  is  shown  in  brackets.  *  p<0.05;  **  p<0.01  versus  controls. 


+/+  x7,800  shi/+  x  7,100  shi/shi  x  18,000 


Figure  4.  Electronmicrographs  of  the  optic  nerves  in 
little,  wild,  and  shiverer  mice.  The  shiverer  mouse  is 
shown  for  comparison  to  a  known  hypomyelinated 
mouse. 

Size  of  white  matter  tracts 

Our  morphometric  analyses  suggest  that  the  little  mice  have  CNS  myelination  of  comparable 
thickness  and  density  to  the  wild-type  mouse.  This  conflicts  with  previously  reported  data  showing 
retarded  neuronal  growth,  an  underdevelopment  of  axons  and  dendrites,  of  layer  5  of  area  6  of 
Caviness  in  the  motor  cortex  (12).  However,  since  the  little  mouse  brain  is  known  to  be  smaller  than 
the  wild  brain,  it  is  possible  that  the  corpus  callosum  and  anterior  commissure  are  smaller  as  well. 
In  order  to  reconcile  this  discrepancy,  the  sizes  of  these  regions  were  compared  by  histological 
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methods.  Comparable  serial  sections  were  compared  in  the  little  and  control  mice  and  the  area  of 
each  region  was  calculated.  One  of  these  sections  in  shown  in  figure  5.  The  area  of  the  anterior 
commissure  of  the  little  mouse  is  calculated  to  be  75.9%  that  of  the  wild  mouse,  and  the  area  of  the 
corpus  callosum  is  82. 1%  that  of  the  wild  mouse. 


Figure  5.  Relative  size  of  the  anterior  commissure  in 
little  mice  and  wild  controls. 


Visual  evoked  potential  testing  of  little  mice 

Although  the  morphology  and  quantity  of  myelin  appears  normal  in  little  mice,  the  function  had  not 
been  determined.  We  selected  a  noninvasive  FVEP  test  procedure  in  order  to  correlate  this  study 
with  the  human  studies  of  18q-  syndrome.  The  averaged  plots  of  the  little  mice  and  wild  controls  are 
shown  in  figure  6.  The  first  positive  peak  represents  the  eye  response  (the  electroretinogram)  which 
was  difficult  to  avoid  recording  since  the  size  of  the  mouse’s  heads  was  so  small.  The  implicit  time 
to  first  negative  peak  and  also  to  the  first  positive  peak  of  the  cortical  response  are  very  reproducible. 

There  is  no  significant  difference  between  the  waveforms  of  the  little  and  control  mice.  We  had 
initially  tested  a  group  of  shiverer  mice  and  their  wild-type  littermates  to  establish  the  feasibility  of 
distinguishing  hypomyelinated  animals  using  this  procedure.  The  first  major  negative  and  positive 
peaks  were  significantly  delayed  in  the  shiverer  mutants  (data  not  shown). 
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Growth  Hormone  Deficient  <G  HD)  Mice  Flash  VEP 


Second* 


Non-Growth  Hormone  D  e  fic  ie  n  t  (NG  H  D  )  M  ic  e  Flash  VEP 


Figure  6.  FEVPs  for  little  (growth  hormone  deficient) 
and  wild  mice. 

Endocrine  characterization  of  shiverer  mice: 

Serum  IGF- 1  and  intrapituitarv  GH  levels  in  shiverer  mice 
No  information  is  available  on  the  impact  of  the  shiverer  mutation  on  the  neuroendocrine  axis. 
Because  of  interest  in  the  potential  effect  of  an  abnormality  of  MBP  on  GH  production,  serum  IGF-1 
of  42-day  old  shiverer  and  control  mice  was  assayed  (Figure  7).  Due  to  the  pulsatile  nature  of  GH 
secretion,  random  serum  GH  determination  does  not  provide  useful  information;  therefore,  IGF-1  was 
chosen  as  an  initial  marker  for  GH  production.  Since  the  number  of  42-day  old  control  animals  is 
small,  statistical  evaluation  could  not  be  undertaken;  however,  a  trend  is  apparent.  Assays  on 
additional  animals  may  make  these  data  significant.  The  functional  significance  of  these  changes  on 
parameters  other  than  somatic  growth  has  not  been  evaluated.  The  results  encouraged  us  to  measure 
intrapituitary  GH  levels,  so  pituitaries  were  harvested  from  60-day  mice  and  assayed  (Figure  7). 
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These  data  indicate  no  significant  difference  in  pituitary  GH  production  and  release  between  the  iittle 


mice  and  normal  controls. 


A 


A 


wild  shi/+  shi/shi  wild  shi/+  shi/shi 


Figure  7.  Serum  IGF-1  levels  (A)  and  intrapituitary 
GH  levels  (B)  in  little  and  control  mice.  Values  are 
means  ±  SD. 


Hybrid  mouse  model: 

Development  of  hybrid  mouse  model 

We  have  demonstrated  that  the  hybrid  {lit  lit ;  -  shi)  is  viable.  EM  analysis  is  being  conducted  to 
establish  the  level  of  myelination  in  these  hybrids  and  to  determine  if  this  differs  from  the  shi 

and  the  lit  lit; — r  mice  during  development.  Unfortunately,  the  level  of  MBP  expression  could  not 
be  evaluated  in  these  hybrids  as  proposed  due  their  strain  differences.  These  mice  are  crosses  of  the 
C57B1  and  C3H  strain.  During  the  course  of  these  studies,  we  discovered  that  the  C3H  strain 
produces  almost  twice  as  much  MBP  as  the  C57B1  strain;  therefore,  crossing  them  results  in  a  large 
amount  of  individual  variability  in  MBP  expression. 

MBP  expression  following  GH  treatment: 

Because  of  the  unexpected  complications  with  measuring  MBP  expression  in  the  hybrid  mouse,  we 
altered  our  experimental  plan.  To  answer  the  question  of  whether  GH  treatment  enhances  MBP 
expression  in  an  animal  that  expresses  less  than  normal  amounts,  we  treated  shi /+  animals  that  express 
only  50%  of  normal  levels  of  MBP.  To  date,  two  groups  of  7  shiverer  heterozygote  mice  have 
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received  daily  GH  injections  for  30  days  commencing  on  postnatal  day  5.  Three  groups  of 
approximately  7  shi/+  mice  have  received  saline  injections  for  the  same  interval.  The  mice  were 
weighed  every  5  days  throughout  the  treatment  phase  and  then  sacrificed  at  the  age  of  65  days. 
Brains  were  weighed  and  harvested  for  RNA  and  protein  analysis.  These  analyses  are  being  initiated. 
Additionally,  one  mouse  from  each  group  was  perfused  and  the  brain  isolated  for  electron 
microscopic  analysis  as  described  for  the  little  mouse  above. 

DISCUSSION 

GH-deficient  little  mice  have  previously  been  reported  to  be  hypomyelinated  (10),  and  this  led  us  to 
propose  the  development  of  a  hybrid  mouse  model  that  is  heterozygous  for  MBP  in  a  growth 
hormone  deficient  background.  We  hypothesized  that  there  may  be  synergistic  effects  between  the 
two  conditions  that  could  result  in  even  less  myelin  production.  To  examine  this,  we  began  by 
•establishing  the  level  of  MBP  expression  and  myelination  in  the  presumed  hypomyelinated  little  mouse 
for  comparison  to  the  hybrid.  We  discovered  that  the  levels  of  two  major  myelin  proteins,  MBP  and 
PLP,  were  essentially  normal  in  these  mutants  during  the  course  of  development.  Moreover, 
morphometric  analysis  of  myelin  sheaths  in  brain  major  white  matter  tracts  showed  no  difference  in 
myelinated  axon  density,  size,  and  myelin  sheath  thickness  between  little  mutants  and  their  wild 
controls.  Histological  analysis  indicated  that  the  size  of  these  white  matter  tracts  is  reduced 
proportional  to  the  reduction  in  brain  size  of  the  little  mouse  as  compared  to  the  wild  controls.  In 
summary,  the  little  mouse  appears  to  have  a  smaller  brain  and  proportionately  smaller  white  matter 
regions  that  are  otherwise  normal  with  regard  to  myelination. 

The  myelin  appears  normal  in  structure  and  quantity  in  little  mice,  but  the  function  remained  in 
question.  Since  auditory  evoked  brainstem  responses  have  been  used  as  a  functional  measure  of 
myelination  in  shiverer  mutants  (27)  as  well  as  humans,  we  initially  attempted  this  test  with  little  mice. 
Unfortunately,  we  discovered  that  ABRs  of  little  mutants  are  abnormal  in  that  signal  processing 
beyond  the  cochlea  is  disrupted.  This  could  be  due  to  structural  abnormalities  of  the  ear  as  a  result 
of  GH  deficiency.  Therefore,  FVEPs  were  conducted  on  the  little  mice  as  a  means  of  measuring  CNS 
function.  Our  results  show  no  functional  difference  as  measured  by  the  FVEP  between  little  mice  and 
controls.  This  further  supports  our  conclusion  of  normal  CNS  myelin  in  the  little  mouse.  These 
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results  were  quite  unexpected  and  provide  interesting  information  regarding  the  impact  of  naturally 
occurring  growth  hormone  deficiency  on  myelination.  A  publication  reporting  these  results  is  in 
progress. 

This  study  also  generated  valuable  results  when  further  characterizing  the  shiverer  mouse.  18q- 
patients  are  haploinsufficient  for  MBP,  but  their  cause  of  GH-insufficiency  is  unknown.  The  GH- 
insufficiency  region  on  chromosome  18  has  been  narrowed  (14),  and  the  MBP  gene  lies  within  this 
region.  To  determine  whether  MBP  could  be  a  candidate  gene  for  GH-deficiency/insufficiency  in 
18q-  syndrome,  we  measured  serum  IGF-1  and  intrapituitary  GH  levels  in  shiverer  mutant  mice  and 
controls.  The  IGF-1  data  indicate  that  homozygous  shiverer  mice  have  reduced  circulating  levels  of 
this  factor  at  42  days.  Since  the  number  of  control  animals  is  small,  statistical  significance  could  not 
be  determined.  In  contrast,  pituitary  GH  levels  in  shiverer  mice  appear  to  be  normal.  Thus,  this  data 
likely  rules  out  MBP  as  a  candidate  gene  for  GH-deficiency  in  18q-  syndrome. 

Development  of  a  hybrid  mouse  model  that  has  only  one  copy  of  the  MBP  gene  in  a  GH-deficient 
background  we  hoped  would  prove  useful  in  examining  the  relationship  between  GH  and  MBP 
production  in  vivo.  We  have  shown  that  this  hybrid  is  viable  and  are  conducting  limited  studies  to 
determine  the  consequences  of  GH-deficiency  on  MBP  haploinsufficiency. 

Treatment  of  heterozygous  shiverer  mice  with  pharmacological  doses  of  GH  during  early  postnatal 
development  may  elucidate  whether  intervention  with  GH  enhances  MBP  production.  It  has  been 
shown  that  GH  has  a  positive  effect  on  myelination,  but  the  specific  effect  on  MBP  expression  is  not 
known.  Two  groups  of  shi/+  mice  have  completed  GH  treatment  and  their  brains  have  been  collected 
for  analysis.  Also,  tissues  from  two  groups  of  saline-treated  shi/+  controls  have  been  harvested.  The 
analyses  are  in  progress. 
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